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Abstract 

We derive a non-local Nambu-Jona-Lasinio model from QCD with a form factor exactly obtained 
in the infrared limit. With this model, having all the parameters properly fixed through those of 
QCD, we consider the case of finite temperature and compute the solution of the gap equation at 
small temperature and small momenta and zero chemical potential. Taking the quark masses to be 
zero, it is proved that the theory undergoes a phase transition with a critical temperature exactly 
determined. These results prove unequivocally that the picture of the vacuum of QCD as a liquid 
of instantons is a very good approximation. 
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I. INTRODUCTION 



Currently, evidences for the existence of a phase transition in QCD, at finite tempera- 
ture and chemical potential, relies on lattice computations. This was firstly realized in a 
pioneering work by Fodor and Katz [l| and further backed up by more recent studies as in 
2 1 and jsl, notwithstanding some criticisms were cast due to the infamous sign problem 
5|. Studies on the behavior of QCD at finite temperature and density, from a theoretical 
standpoint, are generally performed using some phenomenological models like a Nambu- 
Jona-Lasinio model or a sigma model. The appearance of the infamous sign problem in 
attice analysis has prompted some authors to introduce an imaginary chemical potential 
6|-|8| and the need for a consistent agreement between lattice computations and theoretical 
models has prompted the introduction of more general models using the Polyakov loop. In 
this framework some authors were able to prove the existence of a statistical confinement 
of quarks and give account of the phase diagram of QCD j^, [lo|. The relevant point to 
note for our aims is that, when quark masses are taken into account, there is not a real 
phase transition but rather a cross-over between a confined and a deconfined phase and, in a 
same range of temperatures, also a chiral broken symmetry is seen. But when the chemical 
potential and the masses of the quarks are taken to be zero, a first order phase transition is 
indeed expected at a given critical temperature. 

It is clear from this situation that a proof based on first principles of at least a chiral 
breaking of symmetry, starting from the equations of QCD, does not exist yet. Efforts on 



this direction date back to 1980's where chiral perturbation theory come into play |ll|-|l3| 
but did not produce a value for the critical temperature. The difficulty relies essentially on 
our impossibility to obtain a model for the low-energy behavior of QCD directly from the 
theory. Quite recently, we were able to prove that a non-local Nambu-Jona-Lasinio model 



is indeed such a low-energy limit for QCD 



14|-|l6| and this result has been also obtained by 



Kei-Ichi Kondo '17|. The crucial point in our derivation has been an analytical closed form 



for the gluon propagator in the limit of very low energies [18|, Il9 |. 

The form of the gluon propagator is an essential cornerstone result that permits to perform 
explicitly a lot of low-energy computations directly from the equations of QCD. This is 



clearly shown by a recent paper by Hell et al. 



20| . These authors were able to give a complete 



account of a non-local Nambu-Jona-Lasinio model, both at zero and finite temperature, but 
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the form of the propagator was only guessed through the idea that the ground state of QCD 



is that of a hquid of instantons 2l|. We will see below how successful is such a guess as 
our scenario is perfectly consistent with this view. The point is that now, we will see that 
the theory, at zero chemical potential and zero quark masses, indeed undergoes a symmetry 
breaking at low temperature and we are in a position to obtain the critical temperature just 
computed on the lattice. We just note that the value of the critical temperature obtained 
from lattice computation is yet a debated matter as two competing groups obtain not exactly 
the same value. But for our aims it is enough to be in the right range. 

This paper is so structured. In secUl] we present the results on the infrared limit of 
QCD. In sec lIIIW e show again, to make the paper self-contained, a derivation of the non- 
local Nambu-Jona-Lasinio model from QCD. In sec JIVI we solve the gap equation at low 
temperature and low momenta giving the main result. In sec|V]we give the conclusions. 

II. QCD IN THE INFRARED LIMIT 

As usual, our starting point will be the generating functional of QCD. We take 
Sqcd = -\ j d'xW + j d'xj^li^) - ^^y q{x) 

- j d^xic^d^d^c" + ger^'d^A^^'c') (1) 

being F^^ = d^A"^ — d^A^^ + gf°-^'^A^^Al, g the coupling that in this case is dimensionless, 
q{x) are the quark fields, A^(a;) the vector potentials of the Yang-Mills field and the ghost 
field. So, it straightforward to write down 

Our aim is to find a proper approximation in the low-energy limit. We will perform an 
expansion in the inverse of the 't Hooft coupling. In order to manage this functional, it 
appears essential to find a way to reduce this theory to a simpler one. What could make the 
theory manageable is to find a set of classical solutions, in the proper infrared limit of the 
coupling going to infinity, to start a perturbation series for a quantum field theory that holds 
in the same approximation of a strong coupling. A posteriori, we will verify the soundness 
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of our choice of the classical solutions to build up a quantum field theory by comparison 
with numerical solutions on the lattice and for Dyson-Schwinger equations. 



A. Gluon propagator 



With this aim in mind, we have recently proved the following theorem, holding just for 
classical solutions and producing an asymptotic mapping between the scalar field and the 
Yang-Mills theory in the limit of the coupling going to infinity: 



Mapping Theorem. An extremum of the action 

S= i (fx 



is also an extremum of the SU(N) Yang-Mills Lagrangian when one properly chooses A'^^ with 
some components being zero and all others being equal, and A = Ng"^ , being g the coupling 
constant of the Yang-Mills field, when only time dependence is retained. In the most general 
case the following mapping holds 



Al{x)=vy{x)+0{l/VNg) 



being rj'^ constant, that becomes exact for the Lorenz gauge. 



A first proof of this theorem was given in 18| and, after a criticism by Terence Tao, a 



final proof was presented in 



19| also agreed with Tao 22|]. In the following we give a sketchy 



proof for the sake of completeness, but it should be kept in mind that here we are working 
yet with classical solutions . So, let us consider the equation of motion of the scalar field 



0. 



(3) 



Now, we consider a gradient expansion for this equation in the following way. Let us rescale 
the time variable as t — y/Xt. The above equation becomes 



(4) 



and we are in a position to do perturbation theory on this equation in the limit A — ?■ oo 
setting = Yl'^=o -^ "^n- We note at this point a peculiarity of perturbation expansions for 
nonlinear differential equations. Let us consider the small perturbation case and just rescale 
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the field as — )■ ^0 being ^ a function of the coupling A. Applying this rescaling to eq.(|3]) 
we get 

+ AeV' = 0. (5) 
Now we take A' = A^^ and our perturbation expansion is now for A'. So, our weak pertur- 
bation series seems just arbitrary and the same happens to the strong perturbation series 
we are considering, the rescaling in time conspiring to this. Indeed, our strong perturbation 
series is just dual to the weak one and must share the same properties (see |23|). What 
makes inessential this arbitrariness is that the perturbation series must be mathematically 
consistent and one must impose A — )■ for a weak perturbation and A — t- oo for the strong 
one. But these are formally an infinitesimal quantity and an infinite one and multiplying 
them by a constant is inessential: Our expansion will just get an overall multiplying constant 
^. So, at the leading order equation, one has the equation O^^iq + 0q = that admits the 
solution 

00 = /i23sn ^^/it + 6^^ (6) 

being ^ and 9 two integration constants that can be taken depending on space variables. 
But if these are taken exactly constant we have discovered a set of exact solutions of the 
equation we started from. Better, now one can do a Lorentz boost and transform this in a 



covariant set of massive exact solutions 



241 ]. In this case, doing a Lorentz boost boils down 



to a resummation of all the perturbation series in the inverse of A. 

Now, let us consider Yang-Mills equations for a generic SU(N) group and a generic gauge: 

(7) 

As for the scalar field, we implement a strong coupling expansion with the rescaling t — )■ 
\/N gt and imposing the expansion A"^ = Er=o(v^^)'"^?"^- At the leading order of the 
expansion one get the equations 

Now, one can always find a set of components of the Yang-Mills field, chosen to be all equal, 
that reduces this leading order equation to the one of the scalar field 9^^00 + 00 ~ ^ provided 
we do the identification A = Ng^. We can write 

A^^{t,o) = v;<l>{t,o) + o(^-^y (9) 



A Lorentz boost restores covariance and makes r] coefficients depending on momenta. The 
7] coefficients can be fixed through a gauge choice. E.g. in the Landau gauge one can have 

V^vt = Sat (v,. - ^) . (10) 

We do a couple of important considerations. Firstly, we see that gauge invariance is not 
lost due to the mapping and this mapping is an asymptotic one holding in the limit of 
the coupling going to infinity. So, these classical solutions can be used to work in the 
infrared limit in a quantum field theory preserving the substantial physical behavior in the 
ultraviolet limit in the quantum theory for the scalar and the Yang-Mills fields (triviality 
and asymptotic freedom respectively). 

Now, let us evaluate the two-point function for the Yang-Mills field. One gets immediately 



Ki'^-y) = {TA';^{x)Al{y))=v;r]t{mx)<j>{y)) + 0{l/VNg) 

= v;Vu^{x-y) + Oil/VNg) (11) 

where we have set A{x — y) = {T4>{x)(f){y)) for the two-point function of the scalar field. So, 
we need to identify the two-point function for the scalar field in the proper limit. Indeed, 



we have recently proved [25|, |26| that in the limit of the coupling going to infinity, the scalar 
field reaches a trivial infrared fixed point and the two-point function is exactly determined 
as 

n=0 



being 



TT 



2 (_l'jn+l^-in+^)n 



fl~-p"H-i)7p(^ \ . (13) 

being K{i) = , '^^ ^ 1.3111028777, and a formula for the spectrum of the theory, in 

yl+sin^S 

the strong coupling limit, given by 



-Dl^^^^ l - (") 



From the mass spectrum we can identify a string tension that will be useful in the following. 
We set, using the mapping theorem. 
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a = ( ^ ) A = (27rAra,)4A. (15) 
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Here A is an arbitrary parameter arising from the integration of the equations of the theory. 
So, being an integration constant, it should be obtained from experiment. Finally, we note 
;he following functional expansion for the generating functional of the Yang-Mills theory 
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that holds in a strong coupling limit 

A^^ = aJ d'yDf,{x - y)j'"'{y) + 0{f). (16) 

The mapping theorem grants that the propagator entering in this equation is the same given 
in eq. fll2p with a proper choice of the i] parameters. The form of the propagator shared by 



the two theories in the infrared limit, provided a kind of Ka 



len-Lehman representation 



with a non-positive definite spectral function holds (e.g. see (29|), is an evidence that both 
theories are infrared trivial. This can also be seen by the form of the spectrum that has 
only free quasi-p article states but no bounded interacting states. This by no means implies 
that QCD is trivial, rather this theory is infrared safe due to the presence of quarks. 



B. Consistency of the choice of the classical solutions 

So far, we have introduced a set of classical solutions and built up on them a quantum 
field theory without any further support to this choice being the right one or a proof for 
it. Indeed, we have seen that this construction is self-consistent provided we consider a 
increasingly large coupling but we cannot claim that other solutions do exist providing a 
proper description in the same limit or that those gauge configurations grant an optimal 
saddle point for the path integral of the theory. So, the only way we have to be sure that 
this picture is the proper one is to compare it with numerical data. On the lattice, very large 



volumes were considered for the gluon propagator in the Landau gauge in 30l432l| while in 



Ref. 



33 



34j a numerical solution for Dyson- Schwinger equations was provided. 
Our aim will be to show how, increasing the volume, the agreement between numerical 
data and our analytical results tend to coincide. Numerical Dyson-Schwinger equations 
represent our infinite volume limit and we expect a very near coincidence of results in this 
case. 

We consider two kind of lattice computations: A set of volumes till 80^ directly obtained 
with measurements on the lattice for SU(3) and measurements at 128^ recovered from figure 



2 in 



31| for SU(2). We are able to show in this way that, increasing the volume, our 
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propagator describes even more accurately the one measured on the lattice in the deep 
infrared. We would like to point out that the mass gap is different for these two cases as it 
depends on the value of (3 that, just for this section, has nothing to do with temperature but 
is the coupling on the lattice. We see from the figures below that the situation is the one 
depicted with volumes to be taken increasingly large on the lattice to match even better our 
gluon propagator. Note that we consider a weak dependence on the gauge group as showed 



m 



35l ] that is fully consistent with our discussion above. 



p=6 (48)' 

p=6 (64)"^ 

p=6 (80)"^ 

p=6.2 (48)'^ 

p=6.2 (64)"^ 

p=6.2 (80)"^ 
Frasca 




[GeV^] 



FIG. 1: (Color online) Gluon propagator in the Landau gauge for SU(3), 80^ with a mass gap of 
mo = 321 MeV 
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□ Aguilar, Binosi & Papavassiliou 




lO""* 10"^ 10"^ 10"^ 10° 10^ 

[GeV^] 

FIG. 3: (Color online) Gluon propagator in the Landau gauge for SU(3) obtained by numerically 
solving Dyson-Schwinger equations and a mass gap mo = 399 MeV 



This agreement between lattice computations at increasing volume and the perfect match 
for the numerical Dyson-Schwinger equations with our propagator give a strong support to 
our picture and the view that the classical solutions of our choice provide a correct starting 
point for a perturbative quantum field theory in the infrared limit. 



C. QCD in the infrared limit 

Now, we fix the gauge to Landau through eq. lITO!) and change the potential through 
eq. f|T6|) obtaining at the leading order, after noting that the ghost field decouples at this 
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order, 

Sqcd = \j d'xdSnx)D;l{x - y)j'''iy) + J d'xY,q{x) {i$ 



2 



[X) 



- d'y'D;i{x - y')J2i'(y')Y^'i'(y'yj 

q' / 

+Oil/VNg). (17) 
Now, we use the propagator (ITT]) obtaining in the end 
Sqcd = \I d'xdSr''{x)A{x - y)f^{y) + j 

-(7y7^A j d%A{x-y)f^{y)^ q{x) 

-/A f d'xdS'A{x-y')Y,<i{^)^i'<i{^)Y.^'{y')^i,<i'{y') 

q q' 

+0{l/v^g). (18) 



So, we see that the existence of the infrared trivial fixed point in a pure Yang-Mills the- 
ory has the effect to recover, directly from QCD, a non-local Nambu-Jona-Lasinio model 
always reducible to a local one [l^, [l3]- The physics of this model, both at zero and finite 



temperature, has been fully exploited by Hell, Cristoforetti, Roessner and Weise [20|] with 
a substantial difference that they are forced to guess the form of the form factor (the gluon 
propagator) using a model of a liquid of instantons. Here the form factor is directly obtained 
from QCD but we will see below that their guess is excellently good. 



III. NON-LOCAL NAMBU-JONA-LASINIO MODEL AND GAP EQUATION 

In order to define completely the model we have to try to analyze better the propagator. 
We realize from eq.( fT2l) that higher excited states are exponentially damped and so we can 
limit our analysis to a single scalar field interacting with quarks. So, we approximate the 
propagator as A{p) ^ Bq/Ip"^ — rriQ + ie), being mo ~ 1.19^/o' and cr = (0.44 GeV)"^ the string 
tension, and neglect the other contributions coming from higher excited states. This means 
that the Gaussian term | J d'^xd^yj°''^{x)A{x — y)j'^{y) can be rewritten, in this approxi- 
mation using an arbitrary scalar field a that we integrate over, as | J d'^x [{da)^ — mla^], 
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provided we take 



a = v/3(iV2 _ i)/5^A j d^yA{x - y)j{y) 
and writing down the currents as j'^ = •q'^j. So, one finally has 

d'^x q{x) {i$ 



(19) 



Sqcd ~ 2 ^'^^ 



3(iV2 _ 1) 



-/A / d'xdW^ix-y')J2qix)-Yqix)Y,Q'iy')Y^,Q'iy') 

q q' 



+Oil/VNg). 



(20) 



We get a coupling for the a field that can be ignored for our aims. In order to recover in full 
the non-local model of ref. 2o| we have to identify the form factor depending on the gluon 
propagator. We get immediately 



oo 



G 



2" _ (2n + l)2(7r/2K(i))V + ie 2 



C{p) 



(21) 



being G the Nambu-Jona-Lasinio constant that in our case is given by G = 2^(0) = 
^ 0.7854(^Va), so that C(0) = 1, definitely fixed by QCD. 



1^) lla=0 (2n+l)2(7r72if{j))2 



In Ref. 
Ref. 



20|, a guess was put forward for C{p) using a model of liquid of instantons. In 



2l| the form factor for this case takes the form 



Ci{p) 



d 



d'-[io{OMo-hmi{o] 



with 



\p\d 



(22) 



being and Kn Bessel functions. In the following we normalize this function to be 1 at zero 
momenta dividing it by Ci{0). Weise et al. fix the functional form to C{p) = exp{—p'^d'^/2) 
with d~^ ^ 0.56 GeV in order to avoid too much computational weight . We compared 
our C{p) with that given in Ref. 20|, fixing cr = (0.44 GeV)"^ for the string tension and 
rather taking d~^ ~ 0.58 GeV, not much different from Weise et al. guess. The result is 
presented in figJH The agreement is so strikingly good with the instanton form factor that 
our conclusions strongly support a description of the ground state of QCD as an instanton 
liquid. This result was already pointed out in 
for the running coupling in the infrared limit . 



by comparison with lattice results 



36| 
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FIG. 4: (Color online) 
instanton liquid in Ref. 



Comparison between our form factor with that used in Ref. 



2ll | used as a model in Ref. 



2C]. 



201] and the 



With the given expression for the form factor, representing one of the most important 
results given in this paper, we are able to put down the gap equation for massless quarks as 



obtained from ref. 



20|: 

M{p) = C{p)v (23) 



and 

4:NNf f d'^p ^^^^ M{p) 



ml + l/G J {2nY^^^\-' + M^{p) ^^^^ 
being v the v.e.v. of the a field, the number of colors and Nf the number of flavors. Our 
aim is to prove the existence of a phase transition at finite temperature. 
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IV. FINITE TEMPERATURE GAP EQUATION 



We can evaluate the above results at finite temperature by passing to Matsubara sums. 
So, finally we can write down the gap equation as 

M{uk,p) =C{tUk,p)v (25) 

and 



being the Matsubara frequencies Uk = {^n + l)7rT with n an integer. The limits we are 
interested in are those at small momenta and temperature. The first one is needed for 
consistency with the Nambu-Jona-Lasinio model while the second is needed to identify the 
existence of a phase transition. Now, we are in a position to prove the existence of a critical 
point for which v = and the chiral symmetry is restored. Setting v = into eq. fj26l) . we 
have to solve 

ml + l/G ^^^J (27r)3 ujI+P^ 
A possible study of this equation is through numerical techniques. But taking a look at figJU 
after a simple numerical evaluation, we note that the form factor is about 0.8 for a momentum 
of 260 GeV. Lattice computations as those in estimate the critical temperature at about 
170 MeV, well below the limit where the form factor is approximated by unity. This means 
that, for our aims, the form factor can be reduced to a step function cutting at zero at about 
300 MeV and being unity for lower energies. With this crude approximation we are able to 
get an analytic expression for the critical temperature. Indeed, in this case the integral can 
be exactly evaluated and we get 

vrW . 1 



^2 3 

C ^ 9 



(28) 



that proves, starting directly from QCD, that a critical point indeed exists for which chiral 
symmetry is broken. This formula is in close agreement with the one in a recent work by 



Scoccola and Gomez Dumm 37|. The main difference is that we have all fixed through the 
proper value of the mass gap mg due to the form factor. 

Now, we can get an estimation of A, a parameter otherwise fixed by experiment for 
the Nambu-Jona-Lasinio model, just fixing Tc = 0.17 GeV given by lattice computations. 
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Taking a = (0.44 GeVf, g ^ 3, N = 3 and Nj = 2 we get A = 0.77 GeV, a perfectly 
reasonable value for the Nambu-Jona-Lasinio model. This value decreases by increasing the 
number of flavors. So, from this computation we can conclude that both groups in 4| get 
a perfectly reasonable value for the cross-over temperature, providing this can be maintained 
at zero quark masses and chemical potential. 



V. CONCLUSIONS 



We have successfully shown how, starting from the quantum field theory of QCD, this 
theory reduces, in the low-energy limit, to a non-local Nambu-Jona-Lasinio model having all 
the parameters properly fixed to physical values. Besides, the form of the gluon propagator 
in such a low-energy limit is exactly known. This has implied that, when the computation 
is extended to a finite temperature proof of existence of a phase transition, with 

chemical potential and quark masses set to zero, is given. 

This result should only be considered a starting point for future analysis. The most 
important of this is that we have not been able yet to accomodate the Polyakov loop in this 
approach. Presently, this is an important tool for the understanding of the phase diagram 
of QCD. Besides, it will be interesting to see how the cross-over emerges assuming a quark 
mass different from zero and making the gap equation more involved just introducing a 
chemical potential. This is our working program for the near future. 
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